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Abstract 



We examine the uncertainties involved in the ofF-mass-shell extrapolation 
of the K ^ TTTT decay amplitude with emphasis on those aspects that have 
so far been overlooked or ignored. Among them are initial-state interactions, 
choice of the extrapolated kaon field, and the relation between the asymptotic 
behavior and the zeros of the decay amplitude. In the inelastic region the 
phase of the decay amplitude cannot be determined by strong interaction alone 
and even its asymptotic value cannot be deduced from experiment. More a 
fundamental issue is intrinsic nonuniqueness of off-shell values of hadronic 
matrix elements in general. Though we are hampered with complexity of 
intermediate-energy meson interactions, we attempt to obtain a quantitative 
idea of the uncertainties due to the inelastic region and find that they can be 
much larger than more optimistic views portray. If large uncertainties exist, 
they have unfortunate implications in numerical accuracy of computation of 
the direct CP violation parameter e'. 

PACS number(s): 13.20.Eb, 11.55.Fv, 13.75.Lb, ll.lO.Jj 
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I. INTRODUCTION 



Accurate computation of the long-distance QCD corrections for the decay K ^ tctt is 
important in testing the underlying mechanism of CP violation with the e' parameter. When 
one obtains a real amplitude for K ^ mi hj any method of calculation, its validity is limited 
to the unphysical region below vrvr rescattering. On the kaon mass shell one must supply a 
phase by taking account of rescattering However, rescattering affects a magnitude of 

amplitude in general. While the phase of amplitude arises solely from elastic rescattering 
in the case of K ^ irn, the magnitude of amplitude is affected by inelastic scattering as 
well. Many attempts have been made to incorporate the final-state interaction (FSI) in the 
K ^ TTTT decay amplitudes, most recently by Palante and Pich 0] among others. 

The final-state interaction (FSI) theory was first formulated in potential scattering in 
the 1950 's [0, then extended with dispersion relation to relativistic particle physics, often 
using the Omnes-Muskhelishvili (OM) representation |Q. When our interest was in very 
low energy phenomena, it was a good enough approximation to discard all but elastic FSI 
rescattering. In particle physics, however, inelastic rescattering can be potentially important 
to the magnitude of amplitude even when it does not contribute to the phase. If we knew 
the values of the phase at all energies, we could compute the magnitude with the OM 
representation apart from factors of zeros. However, the phase of a decay amplitude above 
the inelastic rescattering threshold has nothing to do with the inelastic scattering phase shift 
of strong interaction. It depends on decay operators even when all quantum numbers are 
the same. 

In addition to inelastic rescattering, strong interactions take place in the initial state 
before K decays. Though it does not contribute to the phase, the initial-state interaction 
(ISI) generates pj^ dependence for the off-shell amplitude just as the FSI does. The existing 
calculations of the FSI largely neglect the ISI effect. Once the inelastic FSI and the ISI are 
included, the OM representation is by no means 'universal" contrary to the statement often 
found in recent literature. 

Extrapolation off mass shell involves a few more basic issues that have not been discussed. 
One is the asymptotic behavior of a decay amplitude far off mass shell, which is directly 
related to the number of zeros and the asymptotic phase S{oo) of the amplitude. In the 
elastic approximation, we may follow potential scattering theory and determine the asymp- 
totic behavior by Levinson's theorem 0. In relativistic physics which contains inelasticity, 
Levinson's theorem does not hold. Therefore we must determine in one way or another the 
value of 6{oo) and the number of zeros and incorporate them in the OM representation. 

Another issue is nonuniqueness of off-shell amplitudes, which is inherent in quantum 
field theory. According to a general theorem |]^, on-shell amplitudes are unique no matter 
what operator one may choose for a particle field as long as it has a correct wave function 
renormalization and a right set of quantum numbers on mass shell. However, off-shell 
amplitudes can depend on choice of particle fields, which is by no means unique in the case 
of hadrons. In particular, the asymptotic behavior is sensitive to the choice. In order to 
write the OM representation uniquely, therefore, we must first define the kaon field and then 
determine the asymptotic behavior and the number of zeros of the decay amplitude. 

Since we think that many of these basic issues in the off-shell extrapolation have been 
either improperly treated or entirely ignored in literature, we attempt to expose them and 
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obtain some quantitative idea of the uncertainties associated with them in this paper. While 
a recent criticism concerns mostly technical aspects of Ref. we focus on more basic 
aspects of the FSI theory 

The paper is organized as follows: In Sec. II, we clarify distinction between the FSI and 
the ISI. For this purpose we write the K ^ Tin amplitude first by the dispersion relation 
for scattering of a weak spurion off an on-shell kaon. In this dispersion relation, the FSI 
generates the s-channel singularities while the ISI contributes to the u-channel singularities. 
Then we write the OM representation by dispersing the kaon mass and show that both 
the FSI and the ISI contribute to the right-hand cuts in the OM representation. In Sec. 
Ill, we relate the number of zeros of amplitude to the asymptotic phase and magnitude 
in the OM representation in general. In Sec. IV, after reviewing the inherent ambiguity 
of an extrapolated particle field in field theory, we study the asymptotic behaviors of the 
K — >■ TTTT amplitudes in QCD by defining the off-shell kaon field by the partially-conserved 
axial-vector current (PCAC) relation. In Sec. V, we interpolate the phases for the dominant 
AI = I amplitude and for the AJ = | amplitude between the inelastic threshold and the 
high energy limit. Then we examine in detail the case that the decay amplitudes have the 
smallest number of zeros including the well-known low-energy zero. In this particular case 
we make quantitative discussion on the extrapolation effect by giving illustrative estimates 
of how much variation occurs in the extrapolation from pj^ = to Tn\. We find that the 
ISI and the inelastic FSI are not negligible even in the case of a single zero. In more general 
cases, uncertainty is larger. The source of the largest uncertainty is in the number of zeros 
and their locations. 

II. DISPERSION RELATIONS 

The decay matrix element for K -kit is written as 

M = {■K{pa)7:{p,r'\HwmK{pK)). (1) 

In actual computation an effective decay operator Ci{fi)Oi is chosen for Hw, which incor- 
porates short-distance QCD effects in by renormalization group. Since the interaction 
above the energy scale /i is included in Cj(/i)Oj, we ought to truncate dispersion integrals 
accordingly, as we shall discuss later. 

A. Spurion scattering 

Treat H^r as the source of spurion S by setting DS" = H]y and write a dispersion relation 
for the on-shell scattering, 

S{p,) + K{pk) 7c{pa) + 7c{p,), {pI = 0) (2) 

with the Mandelstam variables, s = {ps + PkY = iPa + VbYi t = {pk — VaY = iPs — VbYi 
and u = {pk — PbY = {Ps — VaY = ^'k + — s — t. The decay amplitude M(s,t) at 
t < {rriK + '^tt)^ is real analytic in the complex s-plane; M{s*,t) = M{s,t)*. The right- 
hand cut starts at s = 4m^ while the left-hand cut starts at u = {itlk + TnT^)"^, namely, at 
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s = —2771x171^, + ml — t (= sl). (See Fig. 1.) For t fixed small, the s oo limit of M{s,t) 
is determined by Regge asymptotic behavior of K* exchange. Since aK*{fn1) ~ 0.3, one can 
write a once-subtracted dispersion relation in s at the physical point t = as: 

The decay amplitude is obtained by setting = 0, namely, s = mj^ and t = u = m\. The 
first integral in Eq. (j^) due to the s-channel intermediate states describes the FSI, both 
elastic and inelastic. Its discontinuity across the cut is given by 

ImM(s,t)^ = 5]7r(27r)35(p, +p, -p„)(7r|j^(0)|n°"*)(n°"*|i7^(0)|K). (4) 

n 

The lowest two-body inelastic intermediate state is KK. The triangular diagram depicted 
in Fig. 2a is the simplest and probably the most important of the inelastic FSI. In contrast, 
the second integral of Eq. (|^) arises from the u-channel intermediate states which give the 
discontinuity, 

ImM(s,t)i = -^7r(27rf5(p^-p,-pO(7r|i/w^(0)|n°«*)(n°«*|j^(0)|ir). (5) 

n 

The nearest left-hand singularity is generated by the process, K —>■ pK tt-k, namely, 
a virtual dissociation followed by a weak decay. (See Fig. 2b.) The diagram of Fig. 2b 
has no s-channel singularity and therefore cannot be identified with the FSI. Nonetheless, 
it contributes to the s dependence of the amplitude. The process of the K*7i intermediate 
state of Fig. 2c belongs to the same class. We call this class of diagrams as the ISI diagrams. 
Most generally, a single diagram can contribute to both the ISI and the FSI. It is clear from 
the two diagrams in Fig. 2 that the elastic and inelastic FSI alone, namely the s-channel 
singularities, are not sufficient to determine the s dependence of M(s). It is also obvious 
that the phase due to the diagrams of Fig. 2b and 2c are in no way related to scattering 
phase shifts of strong interaction at the corresponding energy. 



B. Phase representation 

Having defined the ISI by the spurion dispersion relation, we now write the OM dispersion 
relation by taking the initial kaon off mass shell. If we follow the reduction formula ||10|| , 
we would reduce one of the pions in addition to the kaon. Rather than handling product 
of three operators j^r, Jk, and H]^, we may take more an intuitive approach by relying on 
perturbative diagram analysis of analyticity [11]. 

Choosing a suitable extrapolation field for the kaon, we can write a dispersion relation 
in the variable pj^ keeping = 0. Since s = {pk + Ps)'^ = Pk for p^ = 0, we denote pj^ by s. 
In terms of the Mandelstam variables previously defined, we fix t and u to keeping s as 
the only variable. Then all singularities appear on the positive real axis of s. (See Fig. 3.) 
The intermediate states KK, pK, and K*7i of the diagrams in Fig. 2 generate cuts starting 
at s = 4m|;, {nip + rriKY , and (m^* + w^tt)^, respectively. To obtain the phase dispersion 
relation of Omnes-Muskhelishvili |^, we write a dispersion relation for the logarithm of the 
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decay amplitude M{s). Keeping in mind that zeros of M{s) are singularities of InM(s), we 
obtain 

where 6{s) = argM(s). In writing the dispersion integral for S{s) with only one subtraction, 
we make the mild assumption that 6{s) does not keep rising or falling indefinitely for large 
s. Though the absorptive parts of the ISI diagrams and the inelastic FSI diagrams both 
appear on the right-hand cuts of the OM representation, they describe physically different 
processes. When we calculate the on-shell decay amplitudes, we are incorporating part of 
the ISI effect through the K ^ tt form factor. To our knowledge, however, nobody has 
included the ISI process in the OM representation. 

III. ZEROS OF DECAY AMPLITUDE 

The large s behavior of M{s) in Eq. is determined by the number of zeros and the 
asymptotic value of the phase 6{oo). It is straightforward to find that, up to logarithmic 
factors, 

lim M(s) ~ s^-'^C'-)/- (7) 

where we have normalized as 5(4m^) = at the threshold. Therefore, once we know the 
large s behavior of |M(s)|, the number of zeros is determined by S{oo). In potential theory, 
M(s) is proportional to the inverse Jost function The Jost function f+{k) |T2| is 

analytic in the upper half plane of the complex variable k{= (s/4 — m^)^/^) and approaches 
unity at \k\ = oo. The zeros of f+{k) mean bound states. If we repeat the argument leading 
to Eq. d^) for f+{k) and use arg/+(A;) = —6{k) and f+{oo) = 1, Levinson's theorem 
results; Nb + -6{oo) = 0, where Ni, is the number of zeros of f+{k), i.e., that of bound 
states.[| In potential theory, therefore, 6{oo) = in any channel that has no bound state. 
When inelastic scattering is included, however, M{s) is no longer as simple as the inverse of 
the Jost function of potential scattering and therefore Levinson's theorem does not hold. 

It should be emphasized that above the inelastic threshold, the phase S{s) of M{s) has 
nothing to do with the phase shift of vrvr — >■ tttt scattering at the corresponding energy. The 
phase of M{s) is a result of complicated interplay of strong and weak interactions unlike 
the scattering phase shift of strong interaction. (See Appendix A.) Therefore, it is nearly 
hopeless to determine 6{oo) either experimentally or directly from theory. All we can hope 
is to determine N — 6{oo)/7c, not and 6{oo) separately, from the large s behavior of |M(s)| 
using Eq. (|^. 

We know that M{s) has one zero in the low-energy region. In the flavor SU(3) symmetry 
limit of strong interaction with no electromagnetic or quark mass difference correction, 
charge-conjugation oddness of parity- violating H\y requires that K ^ nn he forbidden in 



Levinson's theorem is usually quoted as A^'^, = ((5(0) — 5{oo))/Tr. 
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the / = channel [|T3|]. Therefore, the K tttt decay amphtude vanishes at pj^ = Pa — Pb 
in this hmit. It was shown in the soft-pion analysis of the 1960 's that the A J = | amplitude 
of K ^ irn vanishes like ~ 2pj^ ~ Pa ~ Pb 10] ■ If SU(3) breaking enters only through 
the external meson momenta, therefore, M{s) should vanish at s = m^. This low-energy 
off-shell behavior is very robust for the AI = | amplitude. The AI = | decay through 
pure weak interaction [i.e., no electromagnetic or quark mass difference correction) shows 
the same behavior. However, we do not know how many more zeros M{s) may have outside 
the soft-meson region. To study more about zeros, we next look into the large s behavior of 
M{s). 



IV. EXTRAPOLATED FIELD FOR KAON AND ASYMPTOTIC BEHAVIORS 

A well-known theorem of field theory states that no matter what local operator may 
be used for a particle field, the on-shell S-matrix is unique after correct wave-function renor- 
malization is made. The flip side of this theorem is that off-shell matrix elements depend 
on choice of a particle field. We obtain the same on-shell low-energy tttt scattering whether 
we may use the linear cr-model or the nonlinear cr-model. Within the nonlinear model, there 
are different realizations such as the original version by Schwinger , the Callan-Coleman- 
Wess-Zumino (CCWZ) realization |T^, and so forth |T^. All give Weinberg's vrvr scattering 



lengths [^] correctly. In general, however, their off-shell behaviors are different .0 In the 
low-energy expansion, symmetry and kinematical constraints often mask or obscure differ- 
ences among off-shell amplitudes. Nevertheless, the K —>■ im amplitude would depend on 
choice of the extrapolation field if we included the terms of the order higher than O(p^) in 
the expansion. 

S-matrix theory of the 1960's was an attempt to build particle physics theory only with 
on-shell quantities, avoiding off-shell ambiguities. When we used the OM-representation to 
formulate the FSI, however, we actually stepped over the premise of S-matrix theory. Unlike 
the electromagnetic form factors, to which the OM-representation was most successfully 
applied [0, we must take a hadron (the kaon in our case) off mass shell in order to write 
the OM representation for the K ^ nir amplitude. In the case of the electromagnetic form 
factors, the current operator is defined at the beginning. Similarly we must first determine 
what operator is used to describe the off-shell kaon. Only after we have defined the kaon 
field and have chosen QCD as the underlying fundamental theory, do the number and the 
location of zeros as well as the asymptotic value S{oo) become unique. If one chooses one of 
the nonlinear realizations and works with its effective Lagrangians, the large s behavior of 
M{s) would be far more singular and actually meaningless since their applicability is limited 



^ We can see the situation when we express the pion field cj> of the linear (j-model in terms of the 
pion field tt of the nonlinear c-model of CCWZ: cf) = \f7Q,{T^ — 072- (tt • tt) + • • •). The off-shell pion 
amplitude of the linear cr-model is sum of the off-shell one-pion amplitude, the off-shell three-pion 
amplitude and so forth of the nonlinear a-model. Only when the pion is on mass shell, do the 
amplitudes of two models agree with each other. We encounter the same situation when the kaon 
field is taken off mass shell. 
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to low energies. The amplitude M{s) has no unique physical meaning nor direct connection 
to reality except at s = mj^. The phase theorem assures that only the phase of the decay 
amplitude in the elastic rescattering region is unique and physical. Nonuniqueness of off- 
shell quantities should cancel out in the final result of computation of the physical quantity 
M(m|;), but only in principle. We may phrase that the FSI theory formulated with the OM 
representation is not S-matrix theory but field theory even though it uses the technique of 
S-matrix theory. 

Let us study the asymptotic behavior in QCD at large s, that is, at s = 0(/i^) which we 
assume to be the lower end of the asymptopia. Use of the PCAC relation d\A^ = m?pfp(j) 
for the pseudoscalar mesons is appropriate since we are able to work in perturbative QCD 
while maintaining chiral symmetry off shell. We obtain the large s limit of M{s) from 



V2 ^'f^KJK/ J 



We follow the reasoning of perturbative QCD |20| in the approximation of rriu^di'^ fns) 0. 



Since d\A'^ = imsijiLSR—^qL) (<? = u,d), one power of appears from d\A^. Furthermore, 
one insertion of is needed for the right-handed Sr (or s^) of qZsR — 's^qi to turn into Sl 
(or 'sZ) of the weak current. Therefore M{s) must contain ml at large s.Q 

Next we look into chirality matching of the m, d quarks. We first examine the case that 
Oi = Jj^J'^ consists of left-handed currents alone. The leading quark diagram is shown in 
Fig. 4a. In this case, the factor f-^p^ arises from the formation of each pion (the oval blobs 
in Fig. 4). Together with m^, therefore, we obtain the large s behavior of the form 

MM ^ (9) 

up to a power of (Ins)^^. In contrast, when Oi consists of left and right-handed currents as 
the dominant penguin operator Oq does, the pion is formed from a chirality-matched pair 
of WIr ot OnQL with the matrix element rn^fT,/ {rriu + md). This type of decay dominates at 
low energies, e.g., on the kaon mass shell [^. When we go to high energies, however, the 
absence of in the pion formation makes this penguin amplitude softer in s dependence: 

When we combine Eq. (^ or (p!0|) with Eq. (^, we obtain — 5(oo)/7r = 1 or 0. // 
the low-energy zero at s = is the only zero of M(s), then the asymptotic phase should 
approach or vr, depending on the asymptotic behavior of M{s): 

M^} - W (AT = 1; M{s) ^ const). ^^^^ 



^ Though the SU(3) breaking due to one power of m<j is almost compensated by l/mj^, the 
other power of nis remains as a large SU(3) breaking. When we use the PCAC for the kaon field, 
therefore, the amplitude extrapolated to the large s limit shows a large SU(3) breaking through 
the current quark mass difference. 
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The first line applies to the "tree" operator Oi^2, while the second line is for the dominant 
penguin operator Oq. We emphasize that these asymptotic behaviors are specific to the case 
of the PCAC in QCD. They are not necessarily applicable to effective theories for which we 
do not know M(oo). 

For the AI = | decay, the pure weak decay has the zero at s = m^. (Appendix B) If 
this is the only zero (A^ = 1), we have 52(00) = since M{s) ~ s. On the other hand the 
EM penguin matrix element need not vanish at s = m^. The perturbative QCD argument 
presented above remains valid even after one internal photon line is inserted. Therefore 
M{s) ~ constant for the dominant penguin operator Os- The cases of the smallest number 
of zeros for the AI = | amplitudes are, therefore, 

r (A^ = 1; M{s) ~ s) 
62(00) = I (N = 0; Mis) ~ constant), (12) 
[-TT (V = 0;M(s) ~ s). 

The first line applies to the pure weak decay Oi 2, while the second line is for the EM penguin 

If there are more zeros in the amplitude, 5/ (00) shifts upward from those given in Eqs. 
(p!T|) and (p^. The larger the value of |5/(cxd)| is, the more important the contributions from 
the inelastic FSI and the ISI are, as we shall see below. 

V. QUANTITATIVE DISCUSSION 

We first define the problem to which we address in this paper: Suppose that we compute 
with the effective weak interaction Ci{fi)Oi the decay amplitude M{s) at some small value 
So below the tttc threshold. The value of such M{sq) is real though it incorporates all the 
FSI and the ISI that enter through the dispersion integral of S{s). Then we ask how much 
|M(m|')| is different from M(so). We should keep in mind that the value of M(m^) depends 
on choice of the kaon field in principle. Strictly speaking, therefore, this problem makes sense 
only after we have defined the kaon field. 

We focus on the least ambiguous case that M{s) has the minimum number of zeros, 
namely, only one zero at s = for all amplitudes except for the EM penguin and no zero 
for the EM penguin. The zero at or near s = is robust. In contrast, it is rather an ad 
hoc assumption that there is no other zero. We use the values of 5/ (00) predicted by the 
PCAC for the kaon, namely, 5o(oo) = vr for the dominant gluon penguin operator of left and 
right-handed currents and 52(00) = for the AI = | operators, (cf Eqs. (|TI|) and (|T^).) 

When Hw is one of the effective operators Ci{fi)Oi, it is reasonable to match the OM 
representation to the short-distance QCD by cutting off the dispersion integral at s ~ /i^. 
Since /i^ is the spacelike scale of renormalization group, it would be more appropriate to 
make analytic continuation of Cj(/i) to the timelike region. This would generate an imaginary 
part of 0((ln/i)~^) relative to the real part as the high-energy tail of FSI. We can ignore 
it as a small correction if fi is in the perturbative regime of QCD. It has been a subject 
of discussion how far one can lower fi without losing desired accuracy and how one should 
match Ci{fi) to low-energy matrix element computation. If we could choose as low as /i = 1 
GeV, the issues of this paper would mostly disappear since all we would need is the elastic 
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scattering phase shifts. However, such a low cutoff is hard to justify particularly in the 
7 = channel where the scalar resonances /o exist at 1500 MeV and 1700 MeV. While ~ 
2 GeV is probably safe, it may be possible to lower /i close to 1.5 GeV (~ rric)- Near s = /i^, 
6{s) ought to approach its asymptotic value 6{oo). Otherwise the s dependence of M{s) 
computed with the cutoff dispersion integral would not give a smooth asymptotic behavior 
at s = /i^. 

Below the inelastic threshold, the phase of M(s) is equal to the s-wave phase shift of 
elastic tttt scattering by Watson's theorem 0]. The elastic phase shift of / = is positive 
and grows to large values with energy E^, while that of / = 2 stays negative and small 



in magnitude. Though theoretically the inelastic scattering starts at ^/s = 4m7r, in reality 
the inelasticity in the / = channel appears rather suddenly at ^/s = 2mK where the KK 
channel opens. The 1 = scattering phase shift rapidly rises to vr across y/s = 2mK passing 
through the /o(980) resonance just below the KK threshold. As we have emphasized, 
however, the tttt scattering phase shift above ^/s = 2mK has little to do with the phase 
of the K ^ TiTT decay amplitude: Once an inelastic channel opens, the phase theorem no 
longer relates the two phases since the latter depends on the relative sign and magnitude 
of the K — > KK transition matrix element to the K —>■ Tfir decay amplitude too, as we 
shall see below. Therefore the surge of the rni scattering phase shift across ^/s = 2mK does 
not necessarily means that the phase 6o{s) of M{s) rises abruptly in the same way. In fact, 
when an elastic channel strongly couples to an inelastic channel, a partial wave amplitude 
moves quickly toward the center of the Argand diagram and its scattering phase shift often 
undergoes a rapid variation across the inelastic threshold. It is quite possible that the phase 
6o{s) of M{s) varies more slowly across ^/s = 2mK than the tttt scattering phase shift does. 

With the asymptotic values of (5o(oo) = vr and ^2(00) = (cf Eqs. (|TTD and (0), we have 
drawn schematically the behaviors of the phase 5i{s) of M(s) up to s = 2 GeV) in Fig. 
5a by smoothly interpolating between ^/s = 2mx and /i. 



A. Inelastic FSI and ISI 

The phase above the inelastic threshold comes from the ISI and the inelastic FSI. 
For curiosity, we have actually computed the diagrams that generate the lowest singularities 
in the OM representation (Fig. 2a and 2c). The FSI contribution of the diagram in Fig. 2a 
to the absorptive part ImM(s) of /^^ = is: 

ImM^^(s) = (^^) sl^^^I{s)K{s)ReM{s), (13) 

where k{s) is the ratio of the K KK amplitude of A/ = 1/2 to the K — > tttt amplitude 
of AJ = l/2 at Vi, 9l *Kn/'^^ ~ 1.3 is the K*Ktx coupling, and I{s) is a known logarithmic 
function of -^i. We have chosen \ Y.i\(>ut) {out\) + \in){in\) for the intermediate states in 
order to ensure the absorptive part to be real. For a numerical estimate of the magnitude, 
let us choose some value for ^/s, say, ^/s = 2.hmK{— 1-25 GeV), for which I{s) ~ 0.50. We 
set k(s) = (s — m\)/{s — ml), which is the prediction of the octet dominance for Hw and 
SU(3) symmetry for strong interactions. Then we find 

ImM^^~0.39xReM (atv^ = 2.5mK). (14) 
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At higher y/s, the states such as 7777, Oitt, and pp become important intermediate states. 
While rehable estimate is less easy for those decay modes, the inelastic FSI oi K ^ KK 
niT appears sizable. 

Let us turn to the ISI of the diagram of Fig. 2c which generates the lowest ISI singularity. 
There are two independent A J = | transitions for K* pn: 

a(7r ■ p)(KS) + 6(7r x p) ■ (KtS) + h.c, (15) 

where S = (0, 1) represents the AI = ^ spurion. The first term of Eq. (^5|) cancels out in 
the diagram of Fig. 2c by the isospin structure. The second term gives 

ImM^*-(s) = (^ 3K*Krrgp.n \^ (^^) A(s) J(s)ReM(s) , (16) 

where A(s) is the ratio 6/ReM(s), q is the cm. momentum of K*!!, and J{s) is another 
known function of y/s. In contrast to the case of the inelastic FSI of the KK intermediate 
state above, we know little about the value of A(s). For A(s) ^ 1 we obtain at ^/s = 2.bmK 
{J{s) ~ 0.36) 

ImMo^*^ ^ 0.27 X ReMo (at = 2.5mi^), (17) 

The ISI is not negligible either though Eq. (|l^) gives only an order-of-magnitude estimate 
at best. The intermediate states pK, ujK, and so forth start contributing to the ISI at a 
little above y/s = 2.5mK. 

This exploratory numerical exercise indicates that both the ISI and the inelastic FSI are 
large enough to deserve study. Since computation by individual diagrams involves too many 
uncertainties, we instead attempt to extract general trends and some quantitative results 
with minimal assumptions. 



B. / = channel 

We estimate the variation M(m|') — M{so) for the dominant gluon penguin decay with 
So{s) of Fig. 5a. We choose Sq at the SU(3)-symmetry point and quote the variation 
due to the exponential in the vrvr isospin channel I after separating out the factor (s — m^): 

Mj{s) = {s- ml)Mj{s) 

ni=\Mjimj,)\/Mjiml). (18) 

The result is 

ri.50 (^ = 1.5GeV), 
° - \ 1.57 {p = 2 GeV) ^ ^ 

The subtraction in the dispersion integral suppresses the contribution from the high-energy 
tail at 1.5 GeV < < 2 GeV. The authors of Ref. quoted 7^o = 1.41 ± 0.06 without 
including the contribution above the inelastic threshold. We can estimate the effect of the 
ISI and the inelastic FSI by evaluating the phase integral above y/s = 2mK separately from 
the elastic region y/s < 2mK- The contribution from the ISI and the inelastic FSI is: 
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_ f 1.13 (/i = 1.5 GeV), 
~ \ 1.18 (/i = 2 GeV). 



(20) 



The ISI and the inelastic FSI combined contribute to the variation of Mq{s) by a third to 
one half of the elastic FSI. The sizable contribution above the inelastic threshold means that 
we cannot truncate the integral at = 1 GeV and connect to short-distance calculations. 
This is an unfortunate situation since we cannot obtain the values of So{s) at 2m k < a/s < n 
directly from experiment. For comparison we have made an estimate in the case that So{s) 
follows the behavior of the vrvr phase shift across the KK threshold, (see the broken curve 
in Fig. 5a.) The value of TZq^^^^^^^ in this case is: 

^isi+iFsi _ / 1-17 (/i = 1-5 GeV), , . 

'^^ ~ \ 1.23 (/i = 2 GeV). ^^^^ 

Though the ISI and the inelastic FSI are important, their magnitudes are not highly sensitive 
to the s dependence of 5q{s). As we shall see below, it is the value of 5o(oo), namely, the 
number of zeros of the amplitude that they are more sensitive to. 



C. 1=2 channel 



We consider the case of 52(00) = for the pure weak decay {N = 1) and the EM penguin 
decay (A^ = 0). Because of the large enhancement of the A J = | decay over the AI = | 
decay, the 1 = 2 channel of Ks may receive a sizable AI = | contribution from an enhanced 
AI = I decay followed or preceded by long-distance electromagnetic isospin breaking, e.g., 
77-77 mixing.^ Such a contribution may compete with the pure weak decay and the EM 
penguin decay in the AI = | decay. For the AI = | decay generated by long-distance 
electromagnetic corrections, however, 62(3) of M2(s) is no longer equal to the vrvr scattering 
phase shift of J = 2 nor 7 = even below ^/s = 2mK. 

When we ignore the long-distance electromagnetic correction and use 62(3) of Fig. 5a, 
the numbers corresponding to Eqs. (^) and (|20|) are: 

7^2 = 0.91 (/i= 1.5 GeV), (22) 

and 



j^i^si+iFsi ^ 0.99 (^ = 1.5 GeV). (23) 

Our numbers are in perfect agreement with 7^2 — 0.92 ± 0.02 obtained in However 
there is no factor of (s — m^) for the AI = | EM penguin decay. The contribution from 



^In the papers where Cabibbo and Gell-Mann [^] proved the theorem of the vanishing K vrvr 
ampUtude of AI = ^, they even speculated on the possibility that the AI = 3/2 decays could be 
attributed entirely to the electromagnetic correction to the AI = 1/2 transition. The long-distance 
contributions from it — r], it — rj' , and some intermediate states containing a photon were studied 
without FSI m. 
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above 1 GeV is much smaller in J = 2 than in J = because of 62(00) = 0. For the same 
reason, dependence on the value of fi is very mild. For the EM penguin operators Og^io, 
M(s) ~ s so that the asymptotic phase would be 52(00) = — vr (cf Eq. (|12|)). In this case we 
obtain a stronger suppression. It should be emphasized again that Eq. (|^) does not apply 
to long-distance electromagnetic corrections to the gluon penguin operator decay, since we 
cannot equate 62(3) of M2{s) to a tttt scattering phase shift even in the elastic region. The 
OM representation is practically useless in this case. 

D. More zeros (N > 2) 

When Mq{s) has more than one zero, the asymptotic value (5o(oo) is larger. Consider the 
case that Mq{s) has one more zero (A^ = 2). In this case 5o(oo) = 2tt. The second zero must 
be located on the real axis in order to be consistent with real analyticity M{s*) = M{s)*. The 
zero at s = is obtained in the Taylor expansion to 0{p'^) with SU(3) symmetry. When we 
include O(p^), the location of this zero stays as long as SU(3) symmetry is maintained, but 
another zero emerges. If the second zero is outside the soft-meson region, we cannot trust its 
presence. If it is inside or close to the soft-meson region, it means that the next-order terms 
are important at low energies. If the lowest-order chiral Lagrangians give a good description 
of K decays, the second zero should be far from the soft-meson region and therefore the 
second zero factor (s — S2) would not generate very rapid variations at low energies. Our 
hope is that the dominant penguin operator has no second zero and reaches its asymptopia 
quickly. If the second zero exists in the region of \s\ = 0{16tt^Aq^j^), the variation due 
to the second zero factor (s — S2) could be substantial when we extrapolate from s = 
to m\. Lack of our knowledge of its precise location introduces a large uncertainty in the 
extrapolation. 

If we assume that Sq{s) rises linearly in ^/s to 27r from 1 GeV to 2 GeV (Fig. 5b), the 
ratio of the exponential factors at s = m|- and after removing (s — m^)(s — S2) is: 

7^o = 1.64 (/i = 2 GeV). (24) 

In comparison, if So{s) follows the inelastic behavior of tttt scattering phase shift (the broken 
curve in Fig. 5b), the corresponding value could be: 

7^o = 1.83 {fi = 2 GeV). (25) 

In the case of two zeros, therefore, the ISI and inelastic FSI contribution can be as large as 
or even larger than the elastic FSI contribution. 

If M{s) has three zeros, the second and third zeros can appear off the real axis in a 
complex conjugate pair. After we include the variation due to the factor Y\i{s — si), the 
uncertainties get quickly out of control as increases. The uncertainty due to the number 
of zeros and their locations appears to be by far the largest. 

VI. CONCLUSION 

We have studied with the Omnes-Muskhelishvili representation the off-shell kaon mass 
dependence of the K ^ tttt decay amplitudes for the purpose of exposing the uncertainties 
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that have so far not been seriously studied. Determining the phases above the inelastic 
threshold encounters two difficulties. One is complexity of inelastic rescattering and the 
other is physical nonuniqueness of off-shell amplitudes in general. In this paper we have 
studied the off-shell behavior with QCD as the underlying fundamental theory and with the 
PCAC for the extrapolated kaon field. As it is suspected, the largest uncertainty is how 
many zeros exist in the amplitudes and where they are located. We regret that we cannot be 
optimistic in our conclusion as to accuracy or certainty of the off-shell extrapolation. If we 
work in the spurion scattering formalism where all hadrons are on shell, we do not encounter 
the physical nonuniquness due to off-shellness. However, evaluation of the dispersion relation 
is equally formidable since it requires knowledge of the hadr on- spurion scattering with the 
spurion carrying nonvanishing energy-momentum. In either approach, we are unable to work 
only with physical quantities directly measurable in experiment. If we were able to lower 
the cutoff to /i = 1 GeV, we could do without the inelastic region. However, this does not 
seem possible since ^o('S) is not yet asymptotic at 1 GeV. 

If we have very accurate knowledge of M{s) in the low-energy region, we can suppress 
the large s portion of the integral by writing a phase dispersion integral with more than 
one low-energy subtraction ||^. In this case the low-energy values of M{s) that we feed in 
the dispersion relation would incorporate implicitly the contribution from inelastic region. 
The problem would then reduce to determining the physically unique value of M{mj^) with 
the elastic phase shift and with the off-shell values of |M(s)| which are well defined only 
after field theory is specified for mesons. Though the physical nonuniqueness of the off-shell 
|M(s)| should, in principle, cancel out in M(m|-), it is not clear whether we can obtain, 
in practice, accurate enough low-energy values of M{s) from higher-order chiral Lagrangian 
terms. 

Our final comment is on the implications of our study in B decay. The FSI is highly 
inelastic on the B mass shell and the ISI starts at s = {rriB* + Tn^^Y, right above the B 
mass. Unless short-distance strong interactions completely dominate in the FSI and the ISI, 
computation of magnitudes of decay amplitudes would be subject to very large uncertainties, 
far more than in K decay. Several theoretical arguments have been presented in favor of 
short- distance dominance after factorization. Nevertheless, if we should find from experiment 
an FSI phase much larger than the short-distance prediction in some decay mode, our work 
here would be telling us that it is almost hopeless to compute such a decay amplitude with 
decent accuracy. 
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APPENDIX A: PHASE IN THE PRESENCE OF INELASTICITY 

Many of us are aware that the phase of decay amplitude is not determined by strong 
interaction alone at energies where inelastic rescattering occurs. Its value is an outcome of 
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complicated interplay between strong and weak interactions. Sicne there is some confusion 
on this important point in literature, we clarify it in the simplest possible term here. 

Consider the decay amplitude for A — >■ a+b at energies where scattering of a+b is inelastic 
in general. The partial-wave S-matrix of strong interaction scattering is diagonalized by 
eigenchannels: 

= Snn'e^'^". (Al) 
Making time reversal on (n°"*|if^|y4), we obtain for T-invariant Hw 

= {A\Hw\n'n, 

= e''^"(A|i7vK|n°"*). (A2) 

With {A\Hw\n'"'*) = {n°^*\Hw\A)* , this relation leads us to 

M„(s) = M„(s)e*^"(^), (A3) 

where M„(s) is real. This is the phase theorem referred to as Watson's theorem or, to be 
more precise , the Watson- Fermi- Aizu theorem. At energies where only one eigenchannel 
is open in the final state, the phase of a matrix element for any operator having the same 
quantum numbers as H^r is universal and equal to the scattering phase shift. 

However, the situation is different at energies where inelastic channels are open and a 
final state is not an eigenchannel of scattering. Because of T-invariance of strong interaction, 
the transformation between \ab) and eigenchannels \n) can be chosen as orthogonal: 

|a6-*) = Y: Oa,,n\n°^'), (0:,,„ = O,,,.). (A4) 

n 

Therefore the decay amplitude Mab into the observable final state \ab) takes the form 

Mabis) = Y.OabAsWn{s)e''-^'\ (A5) 

n 

The net phase of Mab{s) results from a superposition of the eigenphase shifts weighted with 
the orthogonal transformation matrix elements Oab,n{s) and the weak decay amplitudes 
Mn{s) into the eigenchannels. 

In contrast, the phase shift 6ab of scattering a + b^a + b in the presence of inelasticity 
is defined with the inelasticity T]ab by (a6°"*|a6*") = rjabe^^^''^ Therefore, 

w'^''^^=Y.Olb,ne'''-. (A6) 

n 

The phase 5ab is a quantity determined by strong interactions alone. Above the inelas- 
tic threshold the phase of the decay amplitude M{s)ab has nothing to do with the strong 
interaction phase shift. 



14 



APPENDIX B: LOW-ENERGY ZEROS 

It is straightforward to extend the theorem by Cabibbo and Gell-Mann []T^ to the A/ = | 
decay of ii' — *■ tttt and show that the AI = | amphtude also vanishes in the SU(3) hmit of 
strong interaction. We reconstruct the extension of the proof and show the robustness of 
the low-energy zero. 

When is the current-current interaction in a symmetric product of two octet V — A 
currents, it forms {85} and {27}. Using the SU(3) tensor notation (1,2,3) for {u,d,s), we 
can express the transformation property of Hw as 

Hwr^Hll + H^l (Bl) 

where the upper and lower pairs of indices are both symmetric under interchange. To remove 
{8s} from the right-hand side, we should subtract the trace portion. But it is not necessary 
here since the theorem was already proven for the octet spurion. The relative sign of H^f 
and is positive. In the case of {83} this relative sign aligns Hyy along \q instead of Ay. 
It is crucial to the proof of the theorem for the {27} spurion as well. 

We parametrize the 0^ ^ 0~ + 0~ decay matrix elements of {27} without derivatives in 
terms of the 3x3 octet meson matrix M, instead of U = exp{iXa4>a/ fn) of SU(3) l x SU{3)r: 

M = {MlMiMl + M\MfMl) + (1 ^ 2) + (3 ^ 1) + (1 ^ 2, 3 ^ 1). (B2) 

The relative sign within each bracket is positive, as pointed out above. Otherwise M. would 
point along a wrong component of {27}. 

The parity- violating part of Hy/ is odd under charge conjugation C when Hw is (approx- 
imately) CP invariant. That is, ~Hll under C for H^. It means that the {27} (and 
also {85}) made of symmetric products of the V and A currents must be a C-odd multiplet 
of SU(3).0 Since M M'^ (T = transposed) under C, all terms in Eq. ( p2D transform into 
each other such that Ai —>■ Namely, Ai can make only a C-even {27}. Therefore, we 
cannot construct an SU(3) covariant matrix element having the correct charge conjugation 
property. 

It is obvious in the chiral Lagrangian approach that one cannot write nonderivative 
Lagrangian terms for {8}l nor {27} l of Hw However, chiral symmetry is not really needed 
to prove the theorem; the original proof by Cabibbo and Gell-Mann never used it. 

The theorem fails once SU(3) breaking is included. The decay spurion arising from Hw 
with the electromagnetic interaction or the quark mass difference is no longer a symmetrized 
product of two octets. The EM penguin contains {10}, and {10} too. 

The thereom is valid even after derivatives are included for the mesons: Md^Md'^M 
generates {p2 ■ Ps), which is equal to {pf — pi — pI)/2 and reduces to a common value 
— p^/2 at the SU(3)-symmetry point, pi = p\ = p\ = p^. Therefore all contracted pairs 
of derivatives are factored out as SU(3) invariants so that the SU(3) symmetry property of 
Ma,Ma^M is identical with that of MMM. 



^ The SU(3) charge parity C is defined for a self-charge-conjugate multiplet {n} by the sign of 
|{n};/, l3,y) ib|{n}; /, — Is; — y) under C, or in the tensor notation, T^j. '.' — > ^T^-'''' . 
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Given this theorem, the off-shell K nn amplitude to 0{p'^) is unique up to a propor- 
tionality constant. To prove it, write the off-shell amplitude as a function of p^, pf, and 
by incorporating the Bose statistics for the final TTaTTb in s-wave: 

M = Aipl+pl) + Bpl + 0{p'). (B3) 

In order for }A to vanish at the SU(3) symmetry point, B must be equal to —2A. Conse- 
quently, Ai — —A{2p\ — pi — pi). This argument applies to both {8} and {27}. 

All we need for the proof is M ^ and the fact that A7 = § of is a C-odd {27}. 

Therefore, even if the extrapolation fields of 0" mesons did not respect chiral symmetry off 
mass shell, the theorem would be valid. It is a misnomer to refer to the zero at s = as 
the "chiral zero" . 
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FIGURES 




FIG. 1. Analyticity of the spurion scattering amplitude off K into tttt. 




(a) _ (b) _ (c) 

FIG. 2. (a) An inelastic FSI diagram, (b) The ISI diagram of the nearest left-hand singularity. 

(c) Another ISI diagram, which produces the lowest two-body ISI cut in the OM representation. 




FIG. 3. The s-plane singularity of the OM representation. 




FIG. 4. The quark diagrams for K ^ tttt through (a) the decay operator which contains 
only left-handed currents and (b) the decay operator which contains both left and right currents. 
L(L) and R{R) denote the left-chiral and the right-chiral quark, respectively, of u, d flavors (their 
anti-particles) . 
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FIG. 5. (a) A schematic drawing of the behaviors of Si{s) for 7 = and 1 = 2 with the 
asymptotic values 5o(oo) = tt and 62(00) = 0, respectively. The broken curve is for the case where 
So{s) behaves like the tttt phase shift across s/s = 2mK- (b) So{s) in the case of (5o(oo) = 27r (iV=2). 
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